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The recent discovery of two heavy fermion materials PuCoGa^^ and 
NpPd5Al2\^ which transform directly from Curie paramagnets into supercon- 
ductors, reveals a new class of superconductor where local moments quench 
directly into a superconducting condensate. A powerful tool in the description 
of heavy fermion metals is the large N expansion, which expands the physics in 
powers of about a solvable limit where particles carry a large number (A^) 
of spin components. As it stands, this method is unable to jointly describe the 
spin quenching and superconductivity which develop in PuCoGa^ and NpPd^Al2. 
Here, we solve this problem with a new class of large N expansion that employs 
the symplectic symmetry of spin to protect the odd time-reversal parity of spin 
and sustain Cooper pairs as well-defined singlets. With this method we show 
that when a lattice of magnetic ions exchange spin with their metallic environ- 
ment in two distinct symmetry channels, they are able to simultaneously satisfy 
both channels by forming a condensate of composite pairs between between lo- 
cal moments and electrons. In the tetragonal crystalline environment relevant 
to PuGoGa^ and NpPd^Al2 the lattice structure selects a natural pair of spin 
exchange channels, giving rise to the prediction of a unique anisotropic paired 
state with g-wave symmetry. This pairing mechanism predicts a large upturn 
in the NMR relaxation rate above Tc, a strong enhancement of Andreev refiec- 
tion in tunneling measurements and an enhanced superconducting transition 
temperature Tc in Pu doped Npi^xPuxPd^Al2. 

Large approximations for interacting electron systems have provided an invaluable tool 
for understanding heavy fermion materials [s, and low dimensional magnetism 6, 7]. yet 
to date, they do not encompass heavy fermion superconductivity. We are motivated to 
return to this unsolved problem by the discovery of two new singlet heavy electron super- 
conductors PuGoGa^ and NpPd^Al2 , which transform from Curie paramagnets into 
superconductors without first developing a Fermi liquid. Unlike other heavy electron super- 
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conductors where Cooper pairing is thought to be driven by antiferromagnetic spin fluctu- 



ations 



10 



these materials with a higher transition temperature do not appear to 
be close to a magnetic instability. Moreover, the condensation entropy lost on formation of 
these superconductor is between a quarter and a third of the unquenched spin entropy i? In 2 
of the Curie paramagnet, indicating that the spin-quenching normally associated with the 
Kondo effect is an integral part of the development of superconductivity. 

The difficulty in developing a large N description of heavy fermion superconductivity 
is that the odd time-reversal parity of the electron spin, S — ^ —S, is not preserved by 
SU (N) spins. The inversion of spins under time-reversal protects singlet superconductivity, 
guaranteeing that an electron paired with its time reversed twin is a singlet. However, if 
we extend the theory so that the number of spin components exceeds two, the resulting 
SU{N) spin operators do not all invert under time reversal, so that time-reversed pairs 
of particles cease to be spin singlets. In this paper we show how the time-inversion of 
spins is restored in a new class of symplectic large expansions which link time-reversal 
symmetry to a symplectic property of the electron spin. Our approach enables us to model 
superconductivity on an equal footing with the Kondo effect and leads us to propose that the 
superconductivity in PuCoGa^ and NpPd^Al2 represents a new kind of lattice-coherence 
associated with the Kondo effect, with several observable consequences. 

We begin by demonstrating the link between time reversal and symplectic symme- 
try. Time reversal is an antiunitary operator 6 which acts on an electron wavefunction 
ipcr{x,t) as 6ip{'x,t) = eip*{'x,—t) where ip* is the complex conjugate of the wavefunction 
and iajj = sgn{a)6a,~i3 is the skew symmetric matrix which interchanges "up" and "down" 
amplitudes. 6 is normally written as a product 6 = eK 1^ with the operator K, which 
conjugates quantum amplitudes. Kip = iIj*K. Now consistency of time reversal requires that 
it commutes with spin rotations, U6 = 6U, or 

UeKU^ = eK, (1) 

where U is the unitary rotation operator. But the conjugate of W is its transpose, (W)* = 
U'^, so K converts W into its transpose KU'' = U^K, and consistency requires that 

UiU^ = i. (2) 

This is the symplectic symmetry of spin: a symmetry that must hold if the spin rotation 
group is to remain consistent with the concept of time reversal. The unusual appearance 



of the transpose U'^ , rather than a Hermitian conjugate W reflects the anti-unitary nature 
of time reversal. Replacing U by an inflnitesimal rotation, U = 1 + ia ■ S , we obtain the 
inversion of spins under time reversal S 9S9^^ = i S'^e^ = —S. 

At first sight, this entire set of reasoning immediately extends to the case of fields with any 
even number = 2A; of components, ipa = (^ipi ip-i ■ ■ ■ ipk i^-kj ■ Unfortunately, for N > 2, 
the SU (N) group does not satisfy condition ([2]) for any choice of the matrix e. To preserve 
a consistent definition of both time reversal and spin rotation for > 2, we have to chose 
the symplectic subgroup SP{N), whose elements are actually defined by the consistency 
condition (j2] ) with iajs = sgn{a)6a-i3- Fifteen years ago. Read and Sachdevjl^ made the 
observation that the symplectic group SP{N) allows spin operators that form singlet pairs. 
Their approach has been extensively applied to frustrated magnetism 1^, 1^, the t-J model 



and most recently, to paired Fermi gases 1 131. llTl . |18 |. From this discussion, symplectic 
spins assume an additional importance as the only consistent way to sustain time-inversion 
symmetry in the large A^ limit. 

The spin operators of the SU{N) group are written 



Tap = Ipiipp - j (3) 

where = is the number of particles that make up the spin. Under time reversal. 

Tap a(3T^i3-a, SU{N) spins have no well defined parity. When ijja is a Fermi operator, we 
can also define a charge conjugation operator C that converts particles into holes ipa ccip-a- 
By taking antisymmetric or symmetric combinations of the SU{N) spins with their time- 
reversed version, we may divide them into two sets 

"magnetic" moments Sap = i'ai'p — o^Pip^L/^ip-a, = (~) +) (4) 

which invert under time reversal but are neutral under charge conjugation, and 

"electric" dipoles Vap = i^ai'p + 'S/^'i/'L/j C) = (+, -) (5) 

which are invariant under time reversal, but change sign under charge conjugation[19]. There 
are D = |A^(A^ + 1) magnetic moments that form the generators of SP{N). For the physical 
case of A^ = 2, there are no dipoles, but as A^ becomes large, the SU (N) group contains 
approximately equal numbers of "moments" and "dipoles" 

If we are to sustain the time-reversal properties of spin in a large A^ expansion, we must use 
Hamiltonians that do not contain the unwanted dipole moment operators. This requirement 



delineates the new approach from earher SP{N) treatments of spin systems. Since the 
magnetic moments form a closed SP{N) algebra, their spin dynamics f = -i[H[S], S] 
will remain closed and decoupled from the magnetic dipoles provided the Hamiltonian H 
contains no dipole moments V that spoil the closure. When the local moments are built out 
of fermions, this closure gives rise to a particularly important gauge symmetry. SU{N) spins 
commute with the number operator rij^ at each site, giving rise to a U{1) gauge invariance 
that features heavily in many analyses of correlated electron physics. However, symplectic 
spins also commute with the fermion pair operator = Y2a '^'^j-a'^ja since it is an SP{N) 
singlet, 

[Saf^, ^j] = [SafS, ^]] = [Sap, U,^] = 0. (6) 

In a lattice, these symmetries apply independently at every spin site j, giving rise to an 
S U(2) local gauge invariance. This symmetry was first identified for spin- 1/2 by Affleck et 



al. 20|, who argued for its central role in defining the neutrality of spin. Symplectic spins 
allow us to extend this gauge symmetry to large A^, provided we exclude dipole moment 
operators from the Hamiltonian. This is more stringent requirement than enforcing global 
SP{N) symmetry, and to delineate it from earlier SP{N) approaches, we refer to it as 
"Symplectic- A^" . 

A key distinction between earlier SP{N) approaches and Symplectic-A^, lies in the 
way spin interactions are decoupled. The "dot product" of symplectic spins Si ' S j — 
^Sai3{i)S/3a{j) has a unique decoupling in terms of both particle- hole and singlet pairs: 

5i- 5^2 = -52^521 + 7/^41^21, (W = ±), (7) 

where = a'ipl^ipl^^ creates a valence bond of spins between sites one and two, while 
A21 = J2a "^LV^io- "resonates" valence bonds between sites. The signature 77^ depends on 
whether the field ip is a. boson {rj^p = +1) or fermion {rjjp = — 1). In the large A^ limit, these 



bond variables acquire expectation values. For N = 2, this kind of decoupling [2l| is one of 
many alternative schemes considered by earlier authors. The standard SP{N) decoupling 
procedure omits the last term and is equivalent to a Hamiltonian containing a mixture of 
spins and dipoles, H = JJ2iji^i ' ~ ' ^i)- 



As an initial examination of the Symplectic- A^ procedure we compared its performance 
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when applied to the frustrated Ji — J2 Heisenberg model on a square lattice [23| 



H — Ji Sy: ■ S'x+^ + J2 •S'x ■ Sy:+^', (8) 

where Ji and J2 are the first and next nearest neighbor couplings. We used a Schwinger 
boson representation of the symplectic spins and compared our results with spin wave theory 



and the SP{N) approach. In classical spin wave theory |23l. |2J], the critical spin Sq at 
which quantum fluctuations melt the Neel order in this lattice is strongly dependent on 
the ratio J2/J1, diverging at J2/J1 = 1/2. In SP{N) theory, where the Heisenberg terms 
are expanded purely in terms of pairing fields, the effects of frustration on the critical spin 
are not felt, and the critical spin is independent of J2/J1, an effect presumeably due to 
the presence of dipole spin terms with a ferromagnetic interaction, which over-stabilize the 
antiferromagnet. We find that the inclusion of valence bond terms in Symplectic-N restores 
the strong dependence of the critical spin on frustration. Further discussion of this point is 



found in the online material 22 1. 

We now turn to the app! 
in PuCoGa^ and NpPd^Al2 



ication of Symplectic-A^ to heavy fermion superconductivity 
where the appearance of an SU{2) gauge symmetry in 



the fermionic symplectic-A^ approach has marked physical consequences. These materials 
contain a lattice of local moments, immersed in a sea of electrons. We assume that at low 
temperatures, the Pu and Np ions in these materials behave as Kramer's doublets immersed 
in an electron sea to form a "Kondo lattice". The exchange of spin with its environment 
involves virtual valence fluctuations into ionic configurations with one more, or one less f- 
electron: /" ^ /"^^ =F where n = 3 and 5 for NpPd^Al2 and PuCoGa^ respectively. In 
complex magnetic ions, the partial-wave symmetry of these two virtual processes are distinct 



25l | and we argue that these two processes giving rise to two independent Kondo screening 
channels, F = 1, 2, leading to the following model 

ko- j 

To develop a solvable mean field theory, we examine the family of models where 

Saf^{j) = fU,(^-aPfj_^f,.^. (10) 

are the N component symplectic representation of the magnetic moment at each site j. 
The physical system corresponds to the limit N = 2. The quantities Ji and J2 describe 



two unequal antiferromagnetic exchange couplings in the two symmetry channels while the 
operators ipl^^^ = ['^'r(k)]o(T cj^.^ create electrons in a partial wave state with symmetry 
r = 1, 2, where the matrices and $2k are set by the crystal field symmetry. 

In the tetragonal crystal symmetry environment of NpPd^Al2 and PuCoGa^, the mag- 
netic moments are surrounded by an approximately cubic cage of Pd and Ga ligand atoms, 
respectively. Electrons that interact with the magnetic moment must scatter in one of three 
tetragonal crystal field doublets, labelled {ry,ry,r6}. The states connect selectively 
with the inplane and out-of-plane ligand atoms, as shown in Fig. 1. The Fg crystal field 
state is aligned along the c-axis, overlapping weakly with the nearby ligand ions. This leads 
us to propose a model in which the two F^ channels dominate the spin exchange with the 
f-electrons. 

The presence of two distinct scattering channels plays a central role in our model. For a 
single magnetic Kondo ion, the strongest spin- screening channel always dominates, forming a 
local Fermi liquid of the corresponding symmetry. At Ji = J2 the two channels are perfectly 
balanced, giving rise to a critical state in which the spin screening fluctuates between the 
channels. Many groups have speculated that in a lattice environment, the spins will attempt 
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to avoid this critical state through the development of superconductivity 
Symplectic enables us to develop the first controlled realization of this conjecture, which 
we apply to the new superconductors PuGoGa^ and NpPd^Al2. 



When we expand the Kondo interaction in Q, we obtain Hi = Ylr=i 2 -^'^U) ^ where 



N 



(11) 



describes the spin exchange at site j in channels F = 1,2 and ipjj. = X^k^k^ ^"^^^ creates 
an electron in a Wannier state of symmetry F at site j. This interaction exhibits the local 
SU{2) gauge symmetry — * cosdfu + sui9af'^_^. The important point here, is that this 
symmetry survives for all even N . Earlier efforts have been made to develop SU{2) gauge 



theories of heavy electron systems 



28( 1 and high temperature superconductors [29|], but were 



not justified in terms a controlled expansion. 

When the Kondo interaction is factorized, it decouples into a Kondo hybridization V and 
pairing field A as follows 



{fiVr + a/ _,Ar) + H.C 



A^ 



\Vr\' + \Ar\^ 
Jv 



(12) 



where we have suppressed the site indices j. The mean-field Hamiltonian defined by this 
decouphng becomes exact in the large N limit. Despite the appearance of "pairing terms" 
in the Kondo interaction, the formation of Kondo singlets in a single channel does not lead 
to superconductivity; an SU{2) gauge transformation on the f-electron can always absorb 
the pairing term A into a redefinition of the f-electron: (/J^Vr + af^^A-p) — * Vq/J^. For a 
one-channel Kondo lattice, the Symplectic-A^ and SU (N) large limits are thus identical. 

For two channels, this is no longer the case. Here, it is convenient to define a Nambu 
spinor for the f-electrons and two corresponding matrix SU{2) order parameters 





f)'^ = . ' ^r, = I , (13) 



which transform identically under an SU{2) gauge transformation gj, fj — > gjfj, Vrj 
gjVrj- The product VljVij forms an SU{2) invariant quantity with off-diagonal component 
\& = {VijA2j — V2jAij) : this quantity preserves the local SU{2) invariance, but it breaks 
the global U{1) gauge invariance associated with physical charge, and plays the role of a 
superconducting order parameter. If we carry out an SU{2) gauge transformation that 
removes the Ay, the composite order parameter = VijA2j. These terms lead to Andreev 
reflection off the screened Kondo impurity, as shown in Fig. lb. 

Physically, we may understand this phenomenon as the consequence of the formation of 
a condensate of composite pairs: pairs of electrons, bound to magnetic impurities. When 
an electron scatters off a magnetic impurity, the quantity \1/ determines the amplitude for 



emitting an Andreev hole, leaving behind a composite pair. Detailed analysis 22(] confirms 



this insight and demonstrates the formation of composite order with expectation value 

{^N-2\i^liMlij)S^imN). oc iV,,A2, - V2,A^,). (14) 

In a single impurity model, this order parameter is forbidden, because electrons can never 
change between scattering channels, but in the lattice, electrons travelling between sites no 
longer conserve the channel index, which permits composite order. Once composite order 
develops, the Kondo singlets resonate between the two screening channels, and the resonance 
energy this gives rise to stabilizes the coherent state. 



This basic mechanism was was first proposed in 28|; our large A^ analysis provides the 



first solvable limit in where this mechanism can be rigorously validated, while at the same 



time incorporating the detailed spin-orbit physics of the crystal field-split screening channels. 
When we evaluate the resonant Andreev scattering from the Feynman diagram shown in 
([T]), we find that the conduction electrons acquire a gap that is proportional to the overlap 
of the two hybridization functions Ae(k) oc Tr[$2k'^ik]; as shown in Fig. 1. One of the 
fascinating features of this gap, is that it contains nodes, despite the absence of any such 
nodes in the underlying hybridization. This is a unique consequence of the the hybridization 
between |r~ ) ~ | ± 3/2) and [F^ ) ~ | =1= 5/2) electron states in the composite pair. Since 
these states differ by m = ±4 units of angular momentum, this gives rise to composite pairs 
with I = 4 units of angular momentum that contains gap nodes. 

We have computed the superconducting transition temperature as a function of the 



ratio J2/J1 using Symplectic-A^[22']. Fig. 2 (b) shows the results of a model calculation for a 
two-dimensional Kondo lattice, assuming uniform expectation values for Vr. It is instructive 
to contrast the phase diagrams of the SU{N) and symplectic large limits. In the former, 
there is a single quantum phase transition that separates the heavy electron Fermi liquids 
formed via a Kondo effect about the strongest channel. In the symplectic treatment, this 
quantum critical point is immersed beneath a superconducting "dome" . The Cooper channel 
in the heavy electron normal state guarantees that the secondary screening channel is always 
marginally relevant in the lattice, as first speculated in [2^. This is, to our knowledge, the 
first controlled mean-field theory in which the phenomenon of "avoided criticality" gives rise 
to superconductivity in a model that can be solved exactly. 

There are three concrete consequences of composite pairing that permit our ideas to be 
compared with experiment: 

1. Crystal fields determine the gap symmetry. When an electron scatters between the 
F = F^ channels, it scatters between the |±3/2) ^ |=i=5/2) states, and in so 
doing picks up / = 4 units of angular momentum via spin-orbit scattering. The 
corresponding order parameter then acquires the symmetry of an orbital state with 
/ = 4, and the corresponding order parameter thus has g (1=4)- wave symmetry, with 
eight nodal surfaces, as shown in Fig. 1. The symmetry of this gap is independent 
of the microscopic details, and purely set by the tetragonal crystalline environment of 
the magnetic ions. 

2. Upturn in the NMR. In the approach to the composite ordering transition, the local 
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moments must correlate between sites, and this manifests itself through the develop- 
ment of an enhancement of the NMR relaxation rate. For those systems with maximal 
Tc, the NMR relaxation rate in the normal state is predicted to contain a term derived 
from the interference between the two screening channels, proportional to the product 
TjTf oc Ji(T)J2(T) of the temperature-renormalized Kondo coupling constants. This 
gives rise to an upturn in the NMR relaxation rat e iJ j^ oc [In^ (T/Tc) + vr^]~^, a result 



30|], but which has yet to be tested 



in accord with recent measurements on PuCoGa^ 
in NpPd^Ak. 

3. Andreev reflection. The main driver for this mechanism of heavy fermion superconduc- 
tivity is Andreev reflection off quenched magnetic moments. Conventional Andreev 
reflection involves the direct transfer of an electron from the probe into the pair con- 
densate of the conduction sea. The conventional BTK theory of Andreev reflection 



3l| predicts that such processes are severely suppressed by the large mis-match be- 
tween the probe and heavy electron group velocities. However, in a heavy electron 
system, an electron can also "cotunnel" into the Kondo lattice - a process well-known 
in magnetic quantum dots, whereby the electron flips a localized spin as it tunnels into 



the material 32] ■ In a composite paired superconductor, these processes will result in 
the direct absorption of the electron into the condensate of composite pairs, giving 
rise to an enhanced contribution to the Andreev tunnel current with a Fano resonant 
structure. 

4. Internal Proximity effect. When the Kondo ions in the superconductor are substi- 
tuted by Kondo ions with a larger coupling constant, the robust nature of the gap 
symmetry will protect the superconductor against pair-breaking, and is expected to 
lead to an internal proximity effect, where the Andreev reflection off the substituted 
impurities enliances the superconducting Tg. This effect requires an overlap between 
the gap functions of the two different ions, but this is guaranteed by crystal symmetry, 
providing the same screening channels are operative for both ions. Based on this line 
of argument, we expect that Pu doping of NpPd^Al2 will lead to an enhancement of 
the superconducting Tc. 

Our work currently leaves open, the question of the link between NpPd^Al2 and 
PuCoGa^ and other actinide and cerium systems of similar tetragonal structure, includ- 



ing PuRhGa^ |33| . CeRhIn?.\34i\. Ce/r/ns [35| and CeCo/n^ |36|. Each of these systems 



develops superconductivity, but magnetic susceptibihty measurements indicate that the f- 
moments are more completely quenched at the superconducting transition. We are tempted 
to suggest that that these systems are examples of composite pairing in the parameter range 
where J2/J1 is smaller, and further away from the maximum transition temperature. Neu- 
tron scattering measurements do indeed shown that the two lowest lying crystal field states 
in the Cerium 115 materials are the and Tj states, respectively |32]. Moreover, strong 
Andreev reflections have recently been observed in the tetragonal Ce 115 heavy electron su- 
perconductors, a phenomenon that could be associated with cotunneling into the composite 
pair condensate [38,]. However, at the same time, superconductivity in these systems clearly 
develops in close proximity to antiferromagnetism, so we can not rule out antiferromagnetic 
spin fluctuations as the predominant pairing mechanism in these cases. Future work is 
needed to examine whether the superconductivity observed in these materials is intimately 
linked to the tetragonal crystal field structure, as would be expected in our model. 

In summary, by making the observation of a close link between time-reversal and symplec- 
tic symmetry, we have proposed a new class of large expansion that permits a controlled 
treatment of the Kondo effect and superconductivity on an equal footing. Using this ap- 
proach, we have proposed a model for the heavy electron materials NpPd^Al2 and PuCoGa^, 
which attributes the simultaneous development of superconductivity and spin-quenching to 
the development of a condensate of composite pairs that give rise to resonant bound-states 
that develop a coherent Andreev component to the electron scattering, with a variety of 
experimental consequences. We note in finishing that the methods used in this paper may 
also be applied to one-band electronic systemsflg], an aspect of interest for future work. 
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FIG. 1: (Top) Showing the angular dependence of of the hybridization function I^I'rkP for the 
two tetragonal crystal field scattering states, and Tj with j = 5/2. The states are linear 
combinations of \rja) = cos 6*1 T 3/2) + sin 6*1 ± 5/2) and jPf cr) = cos6l| ± 5/2) - sin6l| T 3/2). 
Figures display the crystal field configurations with maximum overlap with the in-plane and out-of- 
plane ligand atoms, which occurs for 6 ~ vr/lO. (Centre:) Formation of a composite pair between 
the two channels leads to Andreev scattering of conduction electrons between the two channels, 
generating a superconducting gap with symmetry given by the the trace of the overlap between the 
two hybridization functions. (Bottom:) The transfer of four units of orbital angular momentum 
that accompanies this process leads to a gap with / = 4, "g-wave" symmetry, with four large lobes 
of equal sign in the basal plane, as illustrated. 
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FIG. 2: Phase diagram for two channel Kondo lattice, computed in the symplectic large N limit 
for a tetragonal symmetry Kondo lattice in which spin is exchanged in channels Fi = F^ and 
F2 = Fy . The X-axis co-ordinate is the parametric variable x = 2{J2 / J2) / + J2 / Ji) running from 
X = {) to X = 2, corresponding to J2/J1 running from zero to infinity, as labelled. Temperature is 
measured in units of the maximum Kondo temperature of the two channels Tq = max(TR-i, 
Two Fermi liquids of different symmetry develop in the regions of small, or large J2/J1, separated 
by a common region of composite-pairing, delineated by the gray area of the phase diagram. The 
red-point denotes the location of the single-impurity quantum critical point that develops when 
the two channels are degenerate. In the lattice, this point is avoided through the development of 
composite pairing. 15 
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FIG. 3: Upturn in the NMR relaxation rate created by the co-operative interference of the Kondo 
effect in the two channels at different sites, for the extreme case of maximum Tc, where Ji = J2 
(Blue line), compared with measured NMR relaxation rate in PuCoGa^ 30(] (yellow points). Inset 



shows Feynman diagram used to compute this contribution, where dotted lines describe the f- 
fermions, blue lines describe conduction electrons propagating between sites and the curly lines 
describe the Kondo interaction in the particle-hole (red) and particle-particle (blue)^ channels at 
different sites. Temperature is measured in units of the transition temperature Tc [22.] 
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